Excitons in T-shaped quantum wires 
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We calculate energies, oscillator strengths for radiative recombination, and two-particle wave 
functions for the ground state exciton and around 100 excited states in a T-shaped quantum wire. 
We include the single-particle potential and the Coulomb interaction between the electron and 
hole on an equal footing, and perform exact diagonalisation of the two-particle problem within a 
finite basis set. We calculate spectra for all of the experimentally studied cases of T-shaped wires 
including symmetric and asymmetric GaAs/AUGai-z As and InyGai-yAs/AUGai-a, As structures. 
We study in detail the shape of the wave functions to gain insight into the nature of the various 
states for selected symmetric and asymmetric wires in which laser emission has been experimentally 
observed. We also calculate the binding energy of the ground state exciton and the confinement 
energy of the ID quantum- wire-exciton state with respect to the 2D quantum- well exciton for a wide 
range of structures, varying the well width and the Al molar fraction x. We find that the largest 
binding energy of any wire constructed to date is 16.5 meV. We also notice that in asymmetric 
structures, the confinement energy is enhanced with respect to the symmetric forms with comparable 
parameters but the binding energy of the exciton is then lower than in the symmetric structures. 
For GaAs/ALcGai-a; As wires we obtain an upper limit for the binding energy of around 25 meV in 
a 10 A wide GaAs/AlAs structure which suggests that other materials must be explored in order to 
achieve room temperature applications. There are some indications that InyGai-^As/ALjGai-aAs 
might be a good candidate. 
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I. INTRODUCTION 

Optical properties of electrons and holes confined to 
few dimensions are of interest for optical and electronic 
devices. As the dimensionality of the structure is re- 
duced, the density of states tends to bunch together lead- 
ing to a singularity in the ID case. This effect can be 
very useful for low-threshold laser applications. At the 
same time the excitonic interaction in ID is enhanced 
with respect to that in 3D and 2D structures. Quantum 
confinement leads to an increase in the exciton binding 
energy, Ef,, and the oscillator strength for radiative re- 
combination. Both effects provide possibilities for much 
better performance of optical devices such as semicon- 
ductor lasers. 

The binding energy of a ground-state exciton in an 
ideal 2D quantum well is four times that in the 3D bulk 
semiconductor. For the ideal ID quantum wire E\, di- 
verges. This suggests that Ef, for quasi- ID wires can be 
greatly increased with respect to the 2D limit for very 
thin wires with high potential barriers. 3D and 2D exci- 
tons dissociate at room temperature to form an electron- 
hole plasma. To make them useful for real device applica- 
tions, their binding energy needs to be increased and this 
might be achieved by using ID quantum confinement. 

Technologically it is very difficult to manufacture good 
quality ID quantum wires with confinement in both spa- 



tial directions. They can be obtained from a 2D quantum 
well, fabricated by thin-film growth, by lateral structur- 
ing using lithographic methods. The accuracy of this 
method is, however, limited to some ten nanometers and 
thus the electronic properties of samples constructed in 
this way typically have a strong inhomogeneous broad- 
ening. Fortunately it appears possible to achieve quasi- 
1D particles even without a rigorous confinement in any 
of the spatial directions. This has been realised in so 
called V and T-shaped quantum wires. V-shaped quan- 
tum wires are obtained by self-organised growth in pre- 
patterned materials such as chemically etched V-shaped 
grooves in GaAs substrates. The T-shaped quantum 
wire, first proposed by Chang et al. [Q, forms at the 
intersection of two quantum wells and is obtained by the 
cleaved edge over-growth (CEO) method, a molecular- 
beam epitaxy (MBE) technique. The accuracy of this 
method is extremely high and allows fabrication of very 
thin (less than the Bohr radius of an exciton) wires with 
small thickness fluctuations. These structures are cur- 
rently the subject of intensive research and have been 
realised by several groups 0- ||. 

Experimentalists try to optimise the geometry and the 
materials in order to increase the binding energy of the 
excitons, E^, and the confinement energy, E con for pos- 
sible room temperature applications. Up until now, the 
most popular material studied experimentally has been 
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GaAs/Al^Gai-^As. Increasing the Al molar fraction, x, 
should lead to bigger Eb and E con but, unfortunately, 
for larger x the interfaces get rougher which degrades 
the transport properties. Thus optimised geometries for 
lower values of x become more relevant. 

The confinement energy, E con , is the energy difference 
between the lowest excitonic state in the wire and the 
lowest excitonic state in the 2D quantum well. It can be 
directly measured as the difference between the photo- 
luminescence peaks obtained in a quantum wire (QWR) 
and a quantum well (QW). It is, however, not possible 
to measure the exciton binding energy directly. Its value 
has to be obtained from a combination of experimental 
data and one-particle calculations of electron and hole en- 
ergies in a wire. There has been a disagreement between 
the purely theoretical values JlO| |Q and those obtained 
from a combination of experimental data and theoretical 
calculations. The confinement energies, however, tend 
to agree between experiment and purely theoretical cal- 
culations, suggesting that experiment, using combined 
methods where errors tend to accumulate, usually over- 
estimates the binding energy. 

For the 5-nm scale symmetric GaAs/AlAs, Someya et 
al. H reported the largest confinement energy for exci- 
tons in symmetric wires (Figure |]), E con —38 meV and 
£"b=27±3 meV. The largest confinement energy of any 
structure was reported by Gislason et al. |^,|| for their 
optimised wires. Using asymmetric wells with different 
widths and Al content as in Figure ^|, they obtained an 
exciton confinement energy of 54 meV. Recently there 
has also been the first experimental realisation of T- 
shaped wires using Inj,Gai_j,As/Alo.3Gao.7As 0. The 
highest confinement energy reported for this structure is 
34 meV, which is very close to the GaAs/AlAs result, 
and the quality of the structure can be much higher than 
for the GaAs/AlAs case. 

Laser emission from the lowest exciton state in atomi- 
cally smooth semiconductor quantum wires was first ob- 
served by Wegscheider et al. |J in symmetric, T-shaped 
quantum wires made on the intersection of two 70 A 
GaAs quantum wells surrounded by AlGaAs with the Al 
fraction x = 0.35. Recently the same group obtained ex- 
citonic lasing in a 60 A/140 A asymmetric quantum wire 
with a 7% Al filled Stem well (see Figure | ||). They 
reported an interesting observation of two-mode lasing 
in this structure. Under strong excitation they achieved 
simultaneous lasing from two levels in the quantum wire. 
There is a switching between those two lasing modes as 
the temperature or pumping rate is changed. A simple 
rate equation model || gives very good agreement with 
experimental data, which suggests that we have lasing 
from two different states in the quantum wire. 

All calculations published to date which include the 
Coulomb interaction between the electron and hole have 
only examined the ground state exciton. They have used 
either variational methods JT^]- (lj] or other approxi- 
mations p0| , pl| and were performed only for symmetric 
wires and for very limited cases realised experimentally 



in the early days of T-shaped wire manufacturing. With 
the growing experimental realisation of these structures 
as well as the interesting report of lasing phenomena 
there is a need for accurate two-body calculations, treat- 
ing on an equal footing the single-particle potential and 
the Coulomb interaction, of both the ground and excited 
states in the structure. 

Excited states seem to be very important for the op- 
eration of excitonic lasers ||. Calculations of energies, 
oscillator strengths for radiative recombination (i.e, how 
the various states couple to photons) as well as the full 
wave functions for the whole spectra of interest would 
be very beneficial for understanding the origins of cer- 
tain transitions and effects. This could help in the de- 
sign of lasers with better properties and higher maxi- 
mum temperatures for excitonic lasing. The goal is to 
design excitonic lasers which can operate at room tem- 
perature. Also, performing highly accurate calculations 
of the ground state exciton in QWR and the correspond- 
ing QW enables Eb and E con to be obtained for different 
geometries (both symmetric and asymmetric) for a wide 
range of well widths and Al content, x. Such data are of 
great importance for the optimisation of the structures. 

Our method is based on an exact numerical solution of 
the Schrodinger equation in a certain basis within the ef- 
fective mass approximation. The method is not restricted 
to a given number of excited states and we can calculate 
as many of them as required. For some structures we have 
calculated up to 100 excited states. We perform calcula- 
tions for a very wide range of T-shaped wires. In section 
H of the paper the numerical method is discussed in de- 
tail while in section III we present the results. There we 



first study the spectra and wave functions and present a 
discussion of the nature of the various excited states. Fi- 
nally we discuss E con , Eb and the difference between the 
ground-state exciton energy and the first excited-state 
energy, E2-1, as a function of well width Dx and Al mo- 
lar fraction x for the symmetric and asymmetric quantum 
wires. 



II. THE MODEL 

We use the effective mass approximation with an 
anisotropic hole mass to describe an electron in a conduc- 
tion band and a hole in a valence band in the semiconduc- 
tor structures under consideration. The effective mass of 
the hole depends on the crystallographic direction in the 
plane of the T-shaped structure. The electron and hole 
are in the external potential of the quantum wire formed 
at the T-shaped intersection of the GaAs/ALjGai-^As 
quantum wells. The so called Arm quantum well is grown 
in the 110 crystal direction and intersects with a Stem 
quantum well grown in the 001 direction (see Figures |l| 
and ||). In our model the crystal directions 110, 001, 
and 110 correspond to x, y, and z respectively. We con- 
sider symmetric quantum wires where the Arm and Stem 
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well are both of the same width, i.e, Dx = Dy, and 
are made of GaAs. We also consider asymmetric wires 
where the Stem well is significantly wider but filled with 
Al x Gai_ x As with a low Al content to compensate for the 
reduction in confinement energy. Our method is applica- 
ble to any structure regardless of its shape and materials 
provided the external potential is independent of z 

The value of the band-gap is different for the different 
materials used in the well construction. This gives rise to 
the potential barriers at the interfaces between the GaAs, 
Al^Ga^zAs and InGaAs which take different values for 
electrons and holes. In our model the electron and hole 
are placed in external potentials V e (x,y) and Vh(x,y), 
respectively, and interact via the Coulomb interaction. 
We choose the potential in GaAs to be zero and calcu- 
late all potentials in other materials with respect to this 
level. The external potential is independent of z in all 
cases. Sample geometries considered in this work are 
shown in Figs. |l| and [| Using the above model, after the 
separation of the centre of mass and relative motion in 
the z direction, the system is described by the following 
Hamiltonian: 



H 



V 2 V 2 V 2 V 2 

2m e x °< y ' 2m hx Xh 2m hy Vh 2^ z 



V e {x e ,y e ) + V h (x h ,y h ) 



a/ (x e - x h ) 2 + (y e - yh) 2 + z 2 



where z = z e — Zh and — = — H — — . The wave function 

Mz ™c m,hz 

associated with the centre of mass motion in the z direc- 
tion is a plane wave and this coordinate can be omitted 
from the problem. 



A. Numerical method for calculating quantum wire 
exciton states 

We calculate the ground and excited states in the 
structures of interest by a direct diagonalisation method. 
Due to the complexity of the external potential with its 
limited symmetry and sharp edges, none of the standard 
basis sets seem appropriate. We use the following basis 
set: 



ij}(x e ,y e ,x h ,yh,Ze 

E, kir 
Ci,j,kSin(z— 

i ,3,k 



- Z h 
k'TT 



)x e i(x e ,y e ) X Ux h ,yh), (2.1) 



where Xii x e,Ue)/x l j{xh,yh) are electron/hole single- 
particle wave functions for a T-shaped potential without 
the electron-hole Coulomb interaction. In the z direction 
we introduce hard wall boundary conditions and use a 
standing-wave basis set. 

Our basis set does not obey the so called cusp condi- 
tion JlTt which is satisfied whenever two particles come 
together. The divergence in the potential energy when 
the electron and hole come together must be exactly 



cancelled by an opposite divergence in the kinetic en- 
ergy. The exact wave function must therefore have a 
cusp when the electron and hole are coincident. Using a 
basis in which every basis function obeys the cusp con- 
dition would reduce the size of the basis set required. 
For an isotropic hole mass it would be very easy to sat- 
isfy the cusp condition by multiplying the basis functions 
by the factor e ~A\/(^"^) 2 +(^-^) 2 +^ wllich is j ust the 
hydrogenic wave function. Unfortunately there is no an- 
alytical solution when we introduce the anisotropic hole 
mass. Thus we choose not to satisfy the cusp condition 
and therefore have to use a larger basis set. 

The diagonalisation is performed using a NAG library 
routine. Convergence is usually achieved with a basis set 
containing 20 of each of the single-particle wave func- 
tions and 20 standing waves in the z direction. Thus 
20 x 20 x 20 = 8000 basis functions are needed which 
gives 20 6 matrix elements. Only one quarter of the to- 
tal number needs to be calculated as interchanging k\ 
and &2 leaves the matrix element unchanged while inter- 
changing i\ and j\ with ii and j% gives its complex con- 
jugate. This still leaves a great many matrix elements to 
be calculated. Thus to make the calculations feasible the 
matrix elem ents need to be calculated very rapidly (See 
Section II C . 



B. 



Computational method for calculating the 
single-particle wave functions 



The one-particle (electron and hole) wave functions, 
Xi{x e ,y e ) and Xj( x h,yh) in a T-shaped external poten- 
tial are calculated using the conjugate-gradient minimi- 
sation technique with pre-conditioning of the steepest 
descent vector. A detailed explanation of this method 
can be found in reference J16[ . We specify the exter- 
nal potential on a 2D grid and use periodic boundary 
conditions in the x and y directions so that we are able 
to use Fast Fourier Transform (FFT) methods to calcu- 
late the kinetic energy in Fourier space while the poten- 
tial energy matrix elements are calculated in real space. 
The fast calculation of the energy matrix elements is cru- 
cial as they have to be calculated many times during the 
conjugate-gradient minimisation. The FFT provides very 
fast switching between real and Fourier space and makes 
the algorithm much more efficient, but the use of periodic 
boundary conditions introduces the problem of inter-cell 
interactions in the case of two particle calculations. To 
avoid this problem we place the unit cell in the middle of 
another, larger un it ce ll of infinite potential (see Figure 
U and the Section II C). 

We use plane waves as a basis set for the one-particle 
problem. Using this method we can calculate as many 
as 50 states for the electron and 50 for the hole. Very 
good convergence with respect to the number of plane 
wav es and the size of the unit cell is obtained (see Sec- 
tion IniBj). 
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C. Computational method for calculating the matrix 
elements 



The kinetic and potential energies are diagonal in this 
basis and are obtained from the one-particle calculations. 
Thus only the Coulomb matrix elements need to be cal- 
culated. 



A Coulomb matrix element in the basis set (2.1) is a 
5D integral of the following form: 



dx e dy e dx h dy h dzsin(z^- - ^) 



Xt{Xe,Ve)X h i*{xh,yh)- 



n V ~"' "'*' V(Xe ~ X h ) 2 + d/ e - y h Y + Z 2 
■ i k\TT k\TT \ h / \ 

sm{z — ^-)Xh {x e ,ye)x n {Xh, Vh)- 

This integral must be calculated numerically. Numerical 
integration for so many dimensions is very slow and thus 
is not feasible for the case of 20 6 matrix elements. Thus 
another method has to be introduced. 
The above integral is of the form 



dx e dy e dx h dy h dz 

fe(x e ,ye)fh(xh,yh)q(Xe ~ Xh,Ve ~ Vh,z)f z {z). 

Where 

fe(x e ,y e ) = Xll{Xe-,y e )xl 1 {x e ,y e ) 

fh(xh,Vh) = Xizixh^^Xhixh^h), 
q(x e - x h ,y e -yh,z) 



a/ (x e - x h ) 2 + (y e - y h ) 2 + z 2 



Using the Fourier transform and the convolution the- 
orem it can be shown that the above integral is equal 
to: 

j dz J2 F e (~G x ,-G y )*F h (G x ,G y )*Q(G x ,Gy,z). (2 

Where F e , Fh, Q are the 2D Fourier transforms of the 
function f e with respect to x e and y e , fh with respect to 
Xh and yh and q with respect to x e — Xh and y e — yh, 
respectively. Thus the 5D integral can be reduced to a 
ID integral with respect to the z variable and a 2D sum 
in Fourier space. The F e and Fh Fourier transforms can 
be easily calculated using FFTs in real space after mul- 
tiplication of the corresponding xf t (x e , Ve) by xt*(x e ,y e ) 
for electrons and Xi t (xh,yh) by Xi*(xh,Vh)- 

In order to use FFTs we need to introduce periodic 
boundary conditions in the x and y directions as in the 
one-particle calculations. To eliminate interactions be- 
tween particles in neighbouring cells, we place the unit 
cell in the middle of another, bigger unit cell of infinite 
potential (see Figure ||[). 



The distance between the edges of successive small unit 
cells is exactly the width of the small unit cell, L. We 
cut-off the Coulomb interaction at a distance correspond- 
ing to the size of the small unit cell. We therefor consider 
the following form of Coulomb interaction: 



q(x e - x h ,y e ~ yh,z) = 







y/(xe-x h y + (y c -y h y + Z * 



if x e — Xh < L x and 
Ve ~ Vh < L y 
otherwise. 



Particles interact only when their separations in the x 
and y directions are smaller than L x and L y respectively. 
The separations of particles in neighbouring cells is al- 
ways bigger than the cut-off and thus they do not inter- 
act. Particles in the same unit cell are always separated 
by less than that the cut-off distance due to the infi- 
nite potential outside the small unit cell. Thus we take 
into account all of the physical Coulomb interaction and 
completely eliminate the interactions between images. In 
the numerical implementation the infinite potential is re- 
placed by a large but finite potential. Thus the prob- 
ability of the particle being outside the small unit cell 
is effectively zero and we find that the results do not 
depend on the value of this potential for values greater 
than around three times the potential in the ALrGai-zAs 
region. 

The 2D Fourier transform of the 3D Coulomb interac- 
tion with a cut-off cannot be done analytically. Thus we 
put the Coulomb interaction onto a 2D grid as a function 
of relative coordinates x e — Xh and y e — yh for every z 
value. The unit cell in relative coordinates will go from 
—L x to L X1 and —L y to L y respectively. Then for ev- 
ery value of z a 2D FFT is performed with respect to 
X e Xh and y e — yh and the results stored in the 3D ar- 
ray Q(G Xl G y , z). Since this is the same for every matrix 
element the above calculation needs to be performed only 
once. 

The calculations described by Eqn. ^2] need to be per- 
formed for every matrix element. After F e (G x ,G y ) and 
2Fh(G x , G y ) have been calculated the summation over the 
reciprocal lattice vectors G x and G y for every value of z is 
performed. The remaining ID integral in the z direction 
is done numerically, after interpolation of data points, 
using a routine from the NAG library. The dependence 
of the integrand on z is found to be very smooth and thus 
not many points are required to obtain accurate results. 



III. RESULTS 

We perform the calculations for a series of T-shaped 
structures. We calculate energies, oscillator strengths 
and wave functions for the first 20-100 two-particle states 
for symmetric and asymmetric wires. 

For symmetric wires we consider the structure de- 
noted by W which has been experimentally stud- 
ied by Wegscheider et al. M and consists of 
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GaAs/Al .35Gao.65As 70 A quantum wells. Then, keep- 
ing the rest of parameters constant, we vary the quan- 
tum well width from 10 A to 80 A in steps of 10 A in 
order to examine the width dependence of the various 
properties. We also perform calculations for samples de- 
noted by SI and S2 studied by Someya et al. |3|] made 
of GaAs/Al .3Ga .7As (SI) and GaAs/AlAs (S2) quan- 
tum wells of width around 50 A. For the GaAs/AlAs case 
we again vary the well width from 10 A to 60 A. Then 
we take an intermediate value of the Al molar fraction, 
x = 0.56, and vary the well width from 10 A to 60 A 
in order to examine the dependence on the well width 
as well as Al content. Finally we perform calculations 
for 35 A-scale Ino.17Gao.s3As/Alo.3Gao.7As (denoted by 
N4) as well as for 40 A-scale Ino.09Gao.91As/Alo.3Gao.7As 
(denoted by N2) samples as studied experimentally by 
Akiyama et al. H. 

For asymmetric structures we consider the wire studied 
experimentally by Rubio et al. || which consists of a 60 
A GaAs/Al o.35Gao.65As Arm quantum well and a 140 A 
Alo.07Gao.93As/Alo.35Gao.65As Stem quantum well. We 
vary the width of the Arm quantum well from 50 A to 
100 A. We also perform calculations for the asymmetric 
structure studied by a different group j|,D which consists 
of a 25 A GaAsZAlo.3Gao.7As Arm quantum well and a 
120 A Alo.14Gao.86As/Alo.3Gao.7As Stem quantum well. 

In the first part of this section we present the spectra 
for symmetric and asymmetric quantum wires with the 
positions of 2D exciton, ID continuum (unbound elec- 
tron and hole both in the wire) and lDe/2Dh continuum 
(unbound electron in the wire and hole in the well) states 
as well as pictures of representative wave functions. This 
allows us to discuss the nature of the excited states in 
the structures. In the second part we discuss the trends 
in confinement and binding energy and the separation in 
energy between the ground and the first excited states as 
a function of the well width and Al fraction. 

We use a static dielectric constant e=13.2 and a con- 
duction band offset ratio Q c — AE con d/ AE g of 0.65. For 
the difference in bandgaps on the GaAs/AL E Gai_ ;r As in- 
terface we use the following formula: AE g = 1247 x x 
meV for x < 0.45 and 1247 x x + 1147 x (x - 0.45) 2 
meV for x > 0.45. For the electron mass we use m e — 
0.067mo while for the hole mass rrihx = ni} lz = TOmho] = 
0.69 — 0.71mo and rrihy = wimooi] = 0.38mo (mo is the 
electron rest mass). For the Ino.09Gao.91As/Alo.3Gao.7As 
(Ino.17Gao.83As/Alo.3Gao.7As) we use parameters from 
reference [|[: for the electron m e — 0. 0647(0. 0626)mo, 
for the hole rrihy = mhuooi] — 0. 367(0. 358)mo an d 
m hx = m hz = m h[m = 0.682(0.656)m , A£ ff =464(557) 
meV and the band offset was assumed to be 65% in the 
conduction and 35% in the valence band. 



A. Excited states 



1. Symmetric wires 

In Figure || we show spectra (the oscillator strength 
versus energy) for the first 20 (30 in the case of the 70 
A wire) states for the GaAs/Alo.35Gao.65As structure for 
well widths from 10 A to 80 A. A dashed line shows the 
energy of the ID continuum, a dotted line that of the ID 
electron and 2D hole continuum, while the dotted-dashed 
line shows that of the quantum- well 2D exciton. Because 
our system is finite in the z direction, we obtain only a 
sampling of the continuum states; below the continuum 
edge the states are discrete. 

Note that for the experimentally studied 70 A struc- 
ture, the 2D exciton has a lower energy than the com- 
pletely unbound electron and hole in the wire. The sit- 
uation clearly depends on the well width and the cross- 
ing point is between 60 and 70 A. For well widths of 
60 A or smaller, the ID continuum (lDcon) is lower in 
energy that the 2D exciton (2Dexc) with the difference 
being maximal for a width of around 20 A. For widths 
of 70 A or bigger, the lDcon is higher in energy than 
the 2Dexc with the difference growing for increasing well 
width. This effect might be significant for pumping T- 
shaped-wire lasers. Free electrons and holes are excited 
in the whole area of both wells and thus, when the 2D ex- 
citon has a lower energy than the ID continuum, forma- 
tion of the 2D excitons is energetically favourable. These 
excitons can recombine in a well instead of going to the 
wire and forming a ID exciton. Clearly it is more effi- 
cient to have the ID continuum lower in energy than the 
2D exciton. 

By increasing the well width we obtain more states 
that are lower in energy than the lDcon and 2Dexc be- 
ginning with two (ground and the first excited) for the 10 
A well, three for widths between 20-50 A and four states 
for larger widths. 

We now discuss the behaviour of \ij)\ 2 for the 70 A 
case. The wave functions depend on five spatial coordi- 
nates and thus various cuts in 5D space are presented in 
Figures |^ and ^: a) the electron x e , y e position after av- 
eraging over the hole position, b) the hole Xh , yn position 
after averaging over the electron position, and relative co- 
ordinates after averaging over the centre of mass position 
c) the x e — Xh, y e — Vh relative coordinates for z e — Zh = 
and d) the x e — Xh, z relative coordinates for y e — ijh = 0. 

For the ground state we observe that the electron and 
hole are very well localised in the wire with slightly more 
hole localisation. The relative coordinate plots clearly 
show the bound exciton (Figure |5|(1)). 

The electron in the first excited state is localised in the 
wire while the hole already expands into the Arm well. 
The relative coordinate pictures show that the electron 
and hole are bound and form an exciton with an asym- 
metric shape. The size of the exciton is smallest in the 
x direction (the Stem well direction) and the exciton ex- 
pands more into the y (the Arm well where the hole is 
expanded) and free z directions (Figure || (2)). The os- 
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cillator strength of this state is about 1/3 of that of the 
ground state and the state clearly takes the form of a ID 
exciton with its centre of mass in the T-wire. 

It can be seen from the spectra (Figure |J) that there 
are four states (apart from the ground state) with en- 
ergies smaller than lDcon and 2Dcxc. The nature of 
the 3rd and 5th states is very similar to the 2nd one: 
the centre of mass is in the wire and the electron is still 
well localised in the wire while the hole spreads into the 
wells (into both the Arm and Stem wells for the 3rd state 
while only into the Stem well for the 5th one). The rela- 
tive coordinates show the complex, asymmetric shape of 
this excitonic state and the oscillator strength is again 
around 1/3 of the ground state exciton. 

The 4th state with almost zero oscillator strength cor- 
responds to a ID continuum. The electron and hole are 
both in the wire but the relative coordinate pictures show 
an unbound exciton. Within the first 30 states we have 
3 states of that nature: the 4th, 7th and 15th. The 15th 
state is shown in Figure ||: the electron and hole are con- 
fined in the wire (a, b) and there are 3 nodes in the z 
direction and 1 node in the y direction. The other two 
states look similar and differ only in the number of nodes. 
The energy of the 4th state, which is the lowest lDcon 
state, turns out to be lower than the real lDcon obtained 
from our one-particle calculations. This is due to the fi- 
nite size effects. Our method is very well converged with 
respect to the cell size for the bound state and for the 
unbound ones where at least one of the particles is in the 
well. However, for the unbound continuum ID states, 
the particles are very close in the x, y plane because of 
the very small size of the wire and thus the interaction 
is stronger. Consequently it does not decay as fast in 
the z direction as other states and thus we need a much 
bigger unit cell in the z direction to achieve convergence. 
There are however only three such states within the 30 
we examine and we know their true energies from the 
preceding one particle calculations. 

For further excited states up to the 25th, the electron, 
and thus the centre of mass, is still localised in the wire 
while the hole is taking up more and more energetic states 
in both wells, where energies are quantised due to the 
finite size of the cell. Those states can be divided into 
two groups depending on their relative coordinate nature: 
excitonic-like states similar to the second state (Figure || 
(2)) and ionised states like the 22nd which is represented 
in Figure || (22). The oscillator strength of the second 
group is zero (see Figure Wj. 

The 25th state (Figurelq (25)) is the first state where 
the electron is de-localised in both wells, the relative co- 
ordinates and the large oscillator strength shows that it 
is clearly an excitonic-like state. It appears to be a 2D 
quantum well exciton state scattered on the T shaped 
intersection. Its energy is thus higher than that of a pure 
2Dexc. 

The 27th state is the 2D Arm-quantum-well-exciton 
state. It has higher energy than the ground-state 2D 
exciton because the electron and hole wave functions oc- 



cupy higher energy states than the ground state of the 
well due to the presence of the T intersection. 

The 30th state has a very similar nature to the 27th 
but the exciton expands into the Stem instead of the Arm 
quantum well. 

The 25th, 27th and 30th states all have large os- 
cillator strengths (around 3/4 of that of the ground 
state exciton). It is interesting to note that between 
the ground state and those 2D large-oscillator-strength 
states, there is a group of states with relatively low os- 
cillator strengths. The reason for this is that after the 
ground state, there are states where either the wire-like 
electron is bound to the well-like hole and thus they do 
not overlap enough to give big contribution to the spec- 
trum or they consist of a wire like electron with an un- 
bound hole. 

Those quantum-well-like exciton states that scattered 
on the T-shaped potential (like state 25) appear to 
be quite important for the excitonic lasing because of 
their big oscillator strength. In || the authors reported 
two-mode lasing in an asymmetric wire where the laser 
switches between the ground-state exciton and the other 
state whose energy corresponds to the state from the tail 
of the above mentioned states. 



2. Asymmetric wires 

The asymmetric wire that we study in detail consists 
of a 60 A or 56 A GaAsZAlo.35Gao.65As Arm quantum 
well and a 140 A Alo.07Gao.93As/Alo.35Gao.65As Stem 
quantum well ||. The spectrum for the 60 A Arm case 
is shown in Figure |?]. The nature of the states is very 
similar to the case of the symmetric wire. The first two 
excited states are exciton-like and have an electron con- 
fined in the wire while the hole spreads into the well. All 
excited states up to the 20th have the electron confined in 
the wire. The hole spreads to one or both quantum wells 
taking up more energetic states in the well. The rela- 
tive coordinates show either an exciton-like wave function 
(states with nonzero oscillator strength in the spectra of 
Figure |?]) or the case where a hole is confined in the wire 
but is not bound to the electron (states with zero oscilla- 
tor strength in the spectra). Both groups were discussed 
and shown for the symmetric wire. 

The 21st and the 24th states (large oscillator strengths 
in the Figure |?]) have an electron expanding into the Arm 
well. The electron wave function has a node in the wire 
region. The hole wave function spreads into the Arm well 
and has no node for the 21st state and one node in the 
wire region for the 24th state. The relative coordinates 
show the excitonic nature of these states. Thus these 
states correspond to those 2D excitonic states scattered 
on the wire. 

For the asymmetric structure we observe one state 
(the 15th, see Figure ||) which does not correspond to 
any state in the symmetric case. The state is clearly 
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excitonic-like with a large oscillator strength and the rel- 
ative coordinate plots show a very well bound exciton. 
The electron is confined in the wire in the same way as 
the ground state while the hole is clearly ID-like, strongly 
confined in the wire but in a different way. It has a node 
in the wire region. 



B. Trends in confinement and binding energies 



1. Symmetric wires 



We calculate the exciton binding energy, Eb = E e + 
Eh — EiD exc , where E e and Eh are the one-particle en- 
ergies of an electron and a hole, respectively, in the wire. 
We also calculate, using the same method, the exciton en- 
ergy in the quantum well, E 2 Dexc, to obtain the confine- 
ment energy of the ID exciton, E con = E 2 Dexc - Ei Dexc , 
in the wire. 

We perform calculations for a wide range of structural 
parameters. For the GaAs/Al^Gai-^As quantum wire 
we change the well width from 10 A to 80 A for three dif- 
ferent values of the Al content x. The results are shown 
in Figure |J It can be noticed that for a well width bigger 
than 50 A, changing the Al content has very little effect 
on the confinement and binding energies. The difference 
in binding energy between the 60 A GaAs/Alo.35Gao.65As 
and the pure AlAs is only 1.5 meV. Thus it seems more 
promising to change the well width rather than the Al 
content for relatively wide wires. However, for thinner 
wires in the range of 10-50 A, changing the Al content 
is much more profitable then changing the well width. 
The difference in binding energies for 20 A wires with 
Al molar fractions of £=0.3 and x=1.0 is 6.4 meV. This 
increases to 10.6 meV when the width is reduced to 10 

A. 

Eb and E con for Al contents of £=0.35 and £=0.56 both 
approach a maximum for a well width between 10 A and 
20 A. The maximum values for £=0.35 are Ebmax = 17.1 
meV, Econ max — 26.4 meV and for £=0.56 they are 
Ebmax — 19.7 meV, Econ max = 41.4 meV. For the £=1.0 
case, the curve does not have a maximum in the region for 
which calculations has been performed but we consider 
going to wells thinner than 10 A as practically uninter- 
esting. Thus the maximum energies are for Dx=l0 A 
and they are Ebmax — 25.8 meV and Econ m ax = 87.8 
meV. 

E con increases much more rapidly than Eb when the 
well width is progressively reduced. The curves cross for 
a well width between 60 A and 70 A, i.e. for widths of 
60 A or smaller, E con is greater that Eb which means 
that the ID continuum is lower in energy than the 2D 
exciton (as we discussed in section III A 1 ) with the dif- 
ference having a maximum at around 20 A. For widths 
of 70 A or bigger, Eb is greater than E con with the differ- 
ence growing for increasing well width. We also consider 



the difference in energy between the ground state exci- 
ton in the wire and the first excited state as a function 
of the well widths. For the experimentally realised Dx 
— 70 A case, this difference is E%_i = 7.0 meV and the 
maximum value for Dx — 10 A is E 2 -\ m ax — 13.5 meV. 
The maximum value for the GaAs/AlAs at Dx — 20 A 
is 22 meV. 

Although pure AlAs gives the biggest potential offsets 
and thus the biggest binding and confinement energies, 
the GaAs/AlAs interfaces are not very smooth, which 
influences the transport properties. Thus new materi- 
als have to be proposed. Two structures based on In- 
GaAs have been manufactured and measured 0: 35 A- 
scale Ino.17Gao.83As/Alo.3Gao.7As (N4) and 40 A-scale 
Ino.09Gao.91As/Alo.3Gao.7As (N2). The results of calcu- 
lations for these structures are presented in Table || It 
can be seen that energies for the sample N4 are almost 
exactly the same as for the GaAs/AlAs sample S2 sug- 
gesting that these materials might be very good candi- 
dates for structures with large exciton confinement and 
binding energies. 



2. Asymmetric wires 



In order to increase binding and confinement energies, 
the asymmetric T-shaped structure was proposed and re- 
alised by two groups ||,(|,|| ■ 

We calculate E b and E con for the 60 A/ 140 A struc- 
ture with the Stem quantum well filled with 7% Al in 
order to compare with experiment || and then we vary 
the width of the Arm well from 50 to 100 A. One can see 



from Figure 10 that the binding energy is almost inde- 
pendent of the Arm well width in this region, changing 
only from the maximum value of 13.5 meV for Dx = 60 
A to 11.5 meV for Dx = 100 A. The binding energy for 
the 60 A symmetric wire with the same £=0.35 Al mole 
fraction is 13.9 meV - a bit bigger than for the asym- 
metric structure. In contrast, the confinement energy, 
E C on, changes rapidly with the width of the Arm well 
from 4.7 meV for Dx = 100 A up to 33.3 meV for Dx 
= 50 A. For Arm-well widths of 60 A or bigger, the 2D 
quantum-well exciton in the Arm well has a lower en- 
ergy than that for the Stem well, thus the confinement 
energy is calculated with respect to the Arm well exci- 
ton. For the 50 A-wide Arm well, the 2D exciton has 
higher energy than for the Stem quantum well and thus 
the confinement energy is calculated with respect to the 
Stem quantum well. Therefore 33.3 meV is the highest 
confinement energy for this Stem well and changing the 
Arm well would have no effect. Thus the 60 A/ 140 A 
structure is well optimised and its confinement energy, 
E con , is 21.4 meV which is much bigger than that of 14.7 
meV for the 60 A symmetric wire. 

The highest confinement energy so far reported is for 
an asymmetric GaAs/Alo.3sGao.65As wire with a 25 A 
Arm quantum well and a 120 A Stem quantum well filled 
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with 14% Al [||J|] . The experimentally obtained E con for 
this structure is 54 meV. Our calculations however give 
only 36.4 meV which is still the highest among exper- 
imentally obtained structures but much lower than re- 
ported by the authors. Our calculation of E con for five 
different experimentally realised structures agree very 
well with the experimental values and thus it is very 
probable that the value of 54 meV is overestimated. The 
binding energy from our calculations is only 14.6 meV 
for this structure. 

We can conclude from our results that the optimised 
asymmetric structure does not lead to a bigger exciton 
binding energy than the symmetric ones with the same 
parameters. The confinement energy is considerably en- 
hanced and this effect, which can be measured directly, 
has often been used to infer that the binding energy is 
increased. However, our results show that no such re- 
lationship holds between the confinement and binding 
energies. Thus the biggest confinement energy of any 
structure constructed so far of 36.4 meV does not lead to 
the biggest binding energy. Indeed, the binding energy of 
14.6 meV is smaller than the 16.5 meV reported for the 
GaAs/AlAs 50 A-scale symmetric structure || where the 
confinement energy should be only 31.1 meV. It is also 
smaller than expected for a symmetric 25 A-scale struc- 
ture with the same parameters (16.0-16.5 meV). Thus 
asymmetric structures could be useful for applications 
where a large confinement energy is required but appear 
to be less suitable than symmetric wires for applications 
where large binding energies are of interest. 



3. Comparison with experiment and other calculations. 

The comparison between experiment and other pub- 
lished calculations is presented in Table |. The con- 
finement energy of the exciton can be directly measured 
experimentally. Although, due to the strong inhomoge- 
neous broadening of the photolumincsccnce peaks, the 
accuracy of this number is not very high, it is the only 
experimentally proven quantity we can refer to. The ex- 
perimental binding energy needs to be calculated using 
both experimental data and one-particle calculations and 
thus errors might accumulate. Other theoretical methods 
which we refer to obtain the gro und state exciton energy 
using variational techniques fl2|| - JTifl (they differ in the 
form used for the variational wave functions). There are 
also two non- variational calculations for the ground state 
exciton 

Our results for the confinement energy of the ground 
state exciton E con agree very well with experimental val- 
ues for samples SI, N2 and N4 to an accuracy of 1%, 6% 
and 8% respectively. This is indeed very good agreement 
taking into account the strong inhomogeneous broaden- 
ing of the peaks they present. The spectral linewidth of 
the photoluminescence peaks according to the authors is 
around 15 meV which corresponds to a thickness fluctu- 



ation of about 3 A for N2 and N4 @. For the SI and 
S2 samples the authors estimate the experimental error 
due to the inhomogeneous broadening as 2 meV. Agree- 
ment between our calculations and experiment is not as 
good for the S2 sample but for this case additional ef- 
fects are present. For example, AlAs barriers give much 
less smooth interfaces than the lower Al fraction samples 
and this is not taken into account in our model. There is 
also very good agreement (better then 7%) between our 
results and the experimental measurement [^) for asym- 
metric wire R. The earlier E con published by this group 
for the symmetric structure W is probably slightly over- 
estimated. 

There are only two calculations published for the con- 
finement energy. They are based on variational meth- 
ods and were performed only for sample W. Variational 
method 2 [fl3| uses a wave function which takes into ac- 
count correlation in all spatial direction and the agree- 
ment with our results is very good for the confinement 
energy but not so good for the binding energy. 

The variational method proposed by Kiselev et al. |l4| ] 
and denoted here by "3" has a trial wave function which 
has only z dependence in the correlation factor. Their 
binding energy for the sample W differs by only 1 meV 
from our result but their value for the confinement en- 
ergy differs from ours. They perform calculations of the 
binding energy for the whole range of well widths, D x , 
from 10-70 A. This can be compared with our results in 
Figure ^[ Their calculations, like ours, give the maximum 
for Eb and E con for a well width of around 20 A. Their 
binding energy is a bit bigger than the one from our cal- 
culations. They obtained a maximum of Eb = 18.6 meV 
which is 1.5 meV higher than our result. However, their 
confinement energy Econ max = 33.0 meV differs by 7 
meV from our result. Their values of E con are probably 
overestimated. They use the variational technique to cal- 
culate the quantum wire exciton energy but the quantum 
well exciton energy is taken from some other calculations 
of excitons in quantum wells performed using a different 
method and with different parameters, thus errors may 
accumulate. 

The variational method 1 ]l2| ], which uses yet another 
form of trial wave function, has been applied to samples 
SI and S2 to calculate the binding energy Eb. It agree 
quite well with our and other accurate methods. 

The binding energy we obtain shows excellent agree- 
ment with another non-variational calculations by 
Glutsch et al. (IcJ (see Table @). They calculated the 
binding energy only for samples W, SI and S2 and thus 
unfortunately the confinement energy cannot be com- 
pared. The method presented in reference [ pd]| gives much 
lower values for the binding energy than all other meth- 
ods. 

Despite some small differences, all of the theoretical 
methods give much smaller values for Eb than the ex- 
perimental estimates. One has to bear in mind, how- 
ever, that the "experimental" values for Eb (quoted in 
the Table I) are in fact derived from a combination of 
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experimental data and associated theoretical modelling, 
with inherent uncertainties. Our results come from direct 
diagonalisation and are very well converged. Therefore 
we believe that the experimental binding energies are, in 
some cases, considerably overestimated. The real bind- 
ing energy is thus smaller than has been claimed and the 
biggest value for any of the structures manufactured so 
far is 16.5 meV for samples S2 and N4. 



4- Accuracy of the results 

In our method the one-particle energies and wave func- 
tions are calculated first. The one-particle energies are 
very well converged with respect to all the variables such 
as unit cell size, number of points on the grid and the 
number of plane waves to an accuracy of 0.1 meV. We 
use on average as many as 160 000 plane waves which cor- 
responds to 400 x 400 points on the grid (200 x 200 in the 
small unit cell). We obtain excellent agreement between 
our energies for the single electron and hole and those 
obtained by Glutsch et al. PJ. For the 70 A, £=0.35 
symmetric quantum wire we obtain E e — 47.09 meV and 
Eh = 7.47 meV while their results are E e = 47.2 meV 
and Eh = 7.5 meV. According to our calculations there 
is only one electron state confined in the wire and its 
confinement energy E2D-1D (i.e., the difference between 
well-like and wire-like electron states) is 9 meV. This is 
in very good agreement with other methods. L. Pfeifer et 
al. jlfj using eight band k ■ p calculations obtained a con- 
finement energy of 8.5 meV for the same structure. Kise- 
lev et al. juj using the so-called free-relaxation method 
obtained approximately the same value of 9 meV. 

These one-particle wave functions are then used as a 
basis set for the two-particle calculations. The Em exc is 
very well converged with respect to the number of points 
on the grid (as for the one-particle calculations), and with 
the size of the basis set. Convergence is usually achieved 
with about 20 x 20 x 20 (8000) basis functions. In order 
to minimise finite size effects we use quite big unit cells 
(from 43 times the well width, Dx for very thin wires 
(10 A) to 7 times Dx for the 80 A wire). The exciton 
energy Emexc is converged to within about 0.2 meV and 
E-zDexc to within 0.3 meV which gives an accuracy for E\, 
of about 0.3 meV and for E con of about 0.5 meV. 

The other problem which can influence the accuracy of 
the results is the uncertainty associated with the input 
parameters. The electron and hole masses as well as the 
dielectric constant are standard but the potential barriers 
vary a lot depending on the publication. We have found 
quite different values of the potential offsets for the same 
material interfaces in the literature. We have examined 
the influence of this uncertainty on the final results by 
performing calculations for the extrema of the sets of pa- 
rameters found. The binding energy is practically insen- 
sitive to those differences while the confinement energy 
can differ by approximately 2 meV. 



For the parameters that we are using, the results are 
converged to within 0.3 meV for the binding and 0.5 meV 
for the confinement energies. However, one needs to re- 
member that these parameters are not well calibrated 
and this could lead to an additional error in the confine- 
ment energy of about 2 meV. 



IV. SUMMARY 

We have performed an exact diagonalisation within a 
finite basis set of the Hamiltonian which describes an in- 
teracting electron-hole pair in a T-shaped quantum wire. 
We have obtained the ground and excited state energies 
and wave functions for this system. The first group of 
excited states shows an s-like excitonic character where 
the electron is localised in the wire but is bound to the 
hole which spreads into one of the wells. Due to the fact 
that the electron and hole are not localised in the same re- 
gion, we have a group of low oscillator-strength states just 
above the ground state. This group is followed by a num- 
ber of states with large oscillator strength which are 2D 
excitonic states scattered on the T-shaped intersection. 
The excitonic lasing from one of those states has been 
experimentally observed ||. We have also performed a 
detailed study of the exciton binding and confinement en- 
ergies as a function of the well width and Al molar frac- 
tion for symmetric and asymmetric wires. The highest 
binding energy in any structure so far constructed is cal- 
culated to be 16.5 meV which is much smaller than previ- 
ously thought. Our results have shown that for optimised 
asymmetric wires, the confinement energy is enhanced 
but the binding energy is slightly lower with respect to 
those in symmetric wires. For GaAs/Al^Gai-^As wires 
we have obtained an upper limit for the binding energy 
of around 25 meV in a 10 A wide GaAs/AlAs struc- 
ture which suggests that other materials need to be ex- 
plored in order to achieve room temperature applications. 
InyGai^yAs/AlzGai^zAs might be a good candidate. 
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FIG. 2. Shape of the asymmetric T-shaped wire with notations. 
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FIG. 4. Oscillator strength versus energy for the lowest 20 states in a symmetric T-shaped structure for different well widths 
Dx. 
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FIG. 5. Modulus squared of the two-particle wave function for the ground (1), first excited (2) and the 15th (15) state in 
the symmetric T-structure. Electron (a), hole (b) and the relative coordinates x e — Xh, y e — Uh (c), x e — Xh, z (d) probability 
densities are shown. 
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FIG. 6. Modulus squared of the two-particle wave function for the 22th, 25th and the 27th state in the symmetric T-structure. 
Electron (a), hole (b) and the relative coordinates x e — Xh, He — Vh (c), x e — Xh, z (d) probability densities are shown. 
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FIG. 7. Oscillator strength versus energy for the lowest 30 states in an asymmetric T-shaped structure with Dx — 60 A, Dy 
■- 140 A. 
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FIG. 8. Electron (a), hole (b), and the relative coordinates x e — x h , y e — y h (c), x e — x h , z (d) probability densities for the 
15th state in an asymmetric T-shaped structure with Dx — 60 A, Dy = 140 A. 
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FIG. 10. Confinement energy E con — E^Dexc — Emexc, binding energy of the ground state exciton Eb, and the difference 
between the ground and the first excited state E2-1 for an asymmetric wire as a function of the well width Dx, where Dy = 
140 A. 



TABLE I. Binding energy, Eb and the confinement energy fcj CO n — EiDexc — E'2Dexc in meV of the QWR exciton for five 
different samples W, Si , S2 , N2 , N4 obtained from different methods. 
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Var2 j 9.6 11.9 

Var3 k 12 14 



a Sample and experimental values from Ref. 
b Sample and experimental values from Ref. 
c Sample and experimental values from Ref. 
d Sample and experimental values from Ref. 
c Sample and experimental values from Ref. 
Econ is obtained experimentally from the shift between QW and QWR exciton lines. The Eb is obtained indirectly from 
experimental measurement of the QWR exciton line and one-particle calculations. 
s Results of calculations from Ref. 1C|. 
h Results of calculations from Ref. 11]. 
'Results of variational calculations from Ref. J131. 
j Results of variational calculations from Ref. f [L3| . 
k Results of variational calculations from Ref. fll4f . 
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